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y^ , This paper gives a new refined definition of local symmetric order parameters 

QQ ' (OPs)(s-wave, d-wave and p-wave order parameters) of vortex lattice states for sin- 

glet superconductivity.s-wave, d-wave and p-wave OPs at a site (m, n) are defined as 
A, B and E representations of the four fold rotation d at the site (m, n) of nearest 
neighbor OPs {a(m,n)i'^(m+i,n)T) ^tc. The new OPs have a well defined nature such 
that an OP (e.g. d-wave) at the site obtained under translation by a lattice vector 
fj . (of the vortex lattice) from a site (m, n) is expressed by the corresponding OP (e.g. 

d-wave) at the site (m, n) times a phase factor. The winding numbers of s-wave and 
d-wave OPs are given. 
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Growing interest has been focusedjjpn microscopic works on the basis of the Bogoliubov de Gennes 
(BG) equation of an isolated vortextl or a vortex latticea in both conventional and unconventional 
superconductors. In a previous paperEl(referred to as I) we gave a group theoretical classification 
i-j^ ' of vortex lattice solutions of the BG equation in a uniform magnetic field of the extended Hubbard 

r"| . model on a two-dimensional square lattice. In that paper, we considerd conventional definitionsQEliS 

Q ' of s-wave, d-wave and p-wave local symmetric order parameters (OPs), which are valid definitions of 

CJ \ OPs in a spatially uniform superconducting state in the absence of a magnetic field. We showed that, 

e.g., s-wave OP at the site obtained under translation by a lattice vector (of the vortex lattice) from 
a site (m, n) is always a linear combination of s-wave, d-wave and p-wave OPs at the site (to, n). We 
^ I then pointed out that it is not appropriate to assign conventional order parameters at each site. 

VO ' In this paper we give more appropriate definitions of OPs such that an OP at the translated site is 

^£^ , expressed by the corresponding OP at the original site times a phase factor. The model and notation 

used here are the same as those in I. 

We consider singlet superconductivity on a square lattice with inter-atomic spacing a. The Hamil- 
Qv ' tonian for the system is given by 

0\ '. 
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{m,n)ss' 
+ °'{m,n)s'^{™.,n)sO.(^rn,n+l)s''^i™.,n+l)s'\, (1) 

where (m, n) denotes a site in a square crystal lattice, B — (0, 0, B) is a magnetic field, K = 
(7r(/))/(20o)j (ji = Ba? is the magnetic flux through a unit cell of the crystal lattice, and (pQ = ch/2e is 
the flux quantum. The symmetry group of the system is given by 

Go = (e + iC2.)C4rS*, (2) 

where T is the group of the magnetic translation,Ll consisting of the elements T{Maex + NaSy) 
{M,N ^integer ), such that 



and C4 = (e,C^, C2Z, C4],) is the four-fold rotation group around the origin (0,0). Its generator's 
action on the Fermion operator is given by 

Ct-al , =a] , . (4) 

42: [m,n}s { — n,m)s V / 

The action of tC2x (t is time reversal )El on the Fermion operator is 

tC2x ■ \fO-(rn,n)V = ^/ °-{m,-n)i^ 

tC2x-{fa] ,A = f*a] ,,, (5) 

where / is a complex number. Also in Eq. (2), S is the group of the spin rotation (S'[/(2)), and $ is 
the group of the global gauge transformation. 

Hereafter, for illustrative purposes, we restrict our consideration to the symmorphic tetragonal 
vortex lattice states for the case 4> = 4'o/p'^(.P ■ integer), in which the center point of each rotation is 
at a site of the crystal lattice. In this case there are four types of tetragonal vortex lattice solutions 
of the BG equation characterized by invariance groupaa: 

G(i) = {e + tC2x)CLS G = 0,2,±l). (6) 

Here i is a subgroup of T^ consisting elements L{Mpex + NpSy) such that 

L{Mpaex + Npaey)a\^^^^^ = e'i(''^+^'+^^T{Mpaex + Npaey)al^^^^^, (7) 

and 

c'^J^e^^'^Ci, (8) 

i 

where the phase factor e*^'-' should be understood to appear when a group element acts on the space 

VF(a^a^). (The case / = — 1 is overlooked in I.) From L and C invariance of the generalized density 
matrix we have 

{{L{Mpaex + Npaey)a^jn,n}s){L{Mpaex + Npaey)ai^rn'n')s')) 

(e *~ (C4^a(m^„)s)(C4^a(m/„/)s/)) = (a(m,„)sa(m',n')s') (9) 



Then we have 



{0'{m+Mp,n+Np)s O {m'+Mp,n' +Np}s'} 

_ iTT(MN+M+N)+i^(M{n+n')-N{m+m'))i \ 

^ \^{m,n)s^{m,n')s'/i 



\^( — n,m)s^(—7i'.,m')s'/ ^ \^(m.n)s^(m' ,n')s' ) • \^^) 



For the case / — ~l^p ~ 3, using (nOh we obtain a schematic pattern of bond OPs such as 
{0'{m.n)ia{m+i.n)]), which is displayed in Fig. 1. 




FIG. 1. Schematic pattern of bond OP (a(m,n)ifl(m+i,n)T) and (a(m,n)ia(m, n+i)T) for p = 3, I = —1. The 
angle of an arrow with respect to the horizontal x axis represents the phase of the OP. The inner square 
represents a magnetic unit cell. 

It is knowru'El that even in the case of an isolated vortegriine, a d-wave OP induces an s-wave OP. In 
I, we use conventional s-wave, d-wave and p-wave OPsoQ'Q on a site {m, n) as follows: 

S{m,n) = -{{a[,n.n)iO-{rn+l,n)'\) + {o-(m,n)iO-{rn-l.n)^) 
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{m^n)i^{ra,n+l)\ 
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D{m,n) = -{{a(^rn,n)ia{m+l,n)-\) + (a(m,n)ia(m-l,n)T) 
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)-(« 



(m,n)|.^(m,n— l)f 



-i)t)}' 



Px{m,n) — -{(a(m,n)ia(m+l,n)T) ~ (a(m,n)ia(m-l,n)T)}j 
Py{m,n) = ■^{(a(m,n)],a(m,n+l)T) ^ {o-{m,n)i(^{m,n-l)])} ■ 



(11) 



Then we have shown that these OPs at the site {m+Mp, n+Np) are expressed by a linear combination 
of s-wave, d-wave and p-wave OPs at the site {m,n): 

D{m + Mp, n + Np) = }.e'-(MN+M+N)^^f{Mn-Nn^) 

ttN ttM 

X {(cos( — ) - cos( — ))S'(to, n) 

TriV ttM 

+ cos -— + cos(— — £> TO, n) 
2p 2p 

-|-isin(— — )P^(to, n) - isin(— — )P„(TO,n)}, (12) 

2p 2p 

and so on. Then it is very difficult to recognize the symmetry of magnetic translation of the vortex 
lattice by inspecting d wave OPs D{m, n) and so on. Thus it is not appropriate to describe a vortex 
lattice state by assigning the conventional d-wave OP at each site. 

Here we give more refined definitions of s-wave, d-wave and p— wave OPs on a site {m,n). As 
shown in I, 27r/4, tt and (— 27r)/4 rotation around {m^n) (we denote them by C^^{m,n), C2z(to,, n) 
and C^^{m,n)) are given by 



CJ'^{m, n) — T{maex + naey)C'l^T {—maex — naey), 
C2z(to,, n) — T{maex + naey)C2zT {—maex — naey), 
C^^{m,n) — T{maex + naey)C^^T{~maex — naey) 



(13) 



Thus we have 

iy^^[m,n) ■ a(^^j_in)s — e 'a(m,n+l)s: 

4^(m,nj • a(m+i,„)s = e 'a(m.n-i)s, 

= 0,{m.n)s- (14) 

Here we define s-wave, d-wave and p-wave OPs on site(TO, n) as A, B and E representation of the 
four fold rotation group C4{m, n) around the point (m, n) consisting of {e, C^^{m, n), C2z(m, n) and 
C^^{m,n)} as follows : 

S{m,n) = -e*-^"{(a(„,„)|a(„+i,„)|) + {{C2^{m,n)a{^rn.n)i){Ct^{'m,n)a(.^+i^ri)\)) 

+ {{C2zim,n)a(^rn,n)i)iC2z{m,n)a(rn+l,n)^)) 

+ ((C^^(m, n)a(™,„)l)(C^^(m,n)a(„+i,„)|)}}, 
+ ((C2z(m,n)a(„^„)l)(C2^(rn,n)a(„+i,„)|)} 
Px{m,n) = -e'-^"{(a(„,„)|a(„+i,„)|) ~ ((C2z(to, n)a(„,„) J(C2z(m, n)a(„+i^„)-r))}, 

-((C4"^(m,n)a(„,„)^)(C4;(m,n)a(,„+i,„)|)}}. (15) 

Here the factor e*^" is used to obtain a symmetrical expression. Using (llj) we have 

S{m,n) = -{e'- "(a(m^„)^a(m4.i^„)|) + e"' '"(a(m,n)j,a(m,n+i)T) 

I —iKn I \ I iKm / \t 

+ e \a(m,n)Xa(m~l,n)T/ + 6 \a(m,n)J.a()n,n-l)T//j 

D(TO,n) = -{e* "(a(m^„)^a(„i+i^„)|) - e~' '"(a(„^„)j^a(„ „+i)-f) 

+ e ' "(a(m^„)xa(„j_i „)-f) — e* ™(a(m,n)J.Q^()n,n-l)T)}i 
-P:e(?«,«) = ^i*^* "(a(m,ri)ia(m+l,n)T) ^ e"' " (a(m,n)ia(m_l,„)|), 

-e' '"(a(m,n)ia(m,„_l)|)} (16) 

From (|l^) we can see that 

S{m + Mp, n + Np)^ ^^^{MN+M+N)^^fi^In-Nm)g(^^^ ^^^ 
D{m + Mp,n + Np) = e'<MN+M+N)^^^(Mn-Nra)^^^^ ^^ ^ 
K{m + Mp, n + Np)^ ^^^{MN+M+N)^^fiMn-Nm)Jr^^^ „)^ 
Py{m + Mp, n + Np)^ ^^.iMN+M+N)^^fi^In-Nrn)Jr(^^^ „) 

S{-n,m) = e^^''S{m,n), 
D{—n,m) = —e^i^D{m,n), 
P^{-n,m) = -e'^^Py{m,n), 
Py{-n,m) = e^^^Px{m,n) (17) 



Thus these new OPs do not mix due to the invariance magnetic translation L. In Figs. 2 and 3 we 
show schematic patterns of S{rn, n) and D(m, n) corresponding to Fig. 1 for the case I — ~l,p — 3. 
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FIG. 2. A schematic pattern of S(m, n) for the 
case p — 3, I — —1. 



FIG. 3. A schematic pattern of D{m, n) for the 
case p = 3, I = —1. 



TABLE I. Winding number Nw aroud (0,0). 



In the process of making one rotation around the origin (0,0): {m,n) =^ {—n,m) => (— m, —n) => 
(n, — m) => (n, m), total change of the phase of the s-wave OP is 2Trl. The corresponding total change 

of phase for the d-wave OP is 2tt{1 ± 2). Thus 
we obtain winding numbers of s-wave and d-wave 
OPs in Table I. Thus we see that s-wave and 
d-wave have opposite winding numbers for the 
cases I = +1 and I — —1. The opposite winding 
of s-wa«e and d-wave OPs has been pointed by 
VolovikQ in an isolated vortex, which corresponds 
to the case I = —1. Note that this relation holds 
also in the cases of vortex lattice for Z = ±1. 



Nw{S) 









2 


2 


+1 


+1 


-1 


-1 



Nw{D) 

2 



-1 

+ 1 



In conclusion, we obtained well-defined OPs at a site [m,n): S{m,n), D{m,n), Px{m,n) and 
Py{m, n), which have a well-defined nature such that an OP at the site obtained under translation by 
a basis vector (of the vortex lattice) from a site (m, n) is expressed by the corresponding OP at the 
site [m, n) times a phase factor. Similar group theoretical analysis also works for the case in which 
the center of the vortex is located in the middle of a plaquette of the crystal lattice. 
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